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I. INTRODUCTION 



The heavy-heavy-light mesons coupling constants are fundamental objects as they can provide essential information 
on the low energy QCD. Their numerical values obtained in QCD can bring important constraints in constructing the 
meson-meson potentials and strong interactions among them. In the recent years, both theoretical and experimental 
studies on heavy mesons have received considerable attention. In this connection, excited experimental results obtained 
in BABAR, FERMILAB, CLEO, CDF, DO, etc. [I and some physical properties of these mesons have been studied 
using various theoretical models (see for instance 10l-Fl5j ) . 

In this article, we calculate the strong coupling constants, gD*DK* (892) an d gs* bk • (892) m the framework of three- 
point QCD sum rules considering contributions of both D(B) and i^T* (892) mesons as off-shell states, but in the absence 
of radiative corrections. The result of these coupling constants can help us to better analyze the results of existing 
experiments hold at different centers. For instance, consider the B c meson or the newly discovered charmonium states, 
X, Y and Z by BABAR and BELLE collaborations. These states decay to a J/tjj or ip' and a light meson in the final 
state. However, it is supposed that first these states decay into an intermediate two body states containing D q or 
D* with q = u, d or s quarks, then these intermediate states decay into final stats with the exchange of one or more 
virtual mesons. The similar procedure may happen in decays of heavy bottonium. Hence, to get precise information 
about such transitions, we need to have information about the coupling constants between participating particles. 

Calculation of the heavy-heavy-light mesons coupling constants via the fundamental theory of QCD is highly 
valuable. However, such interactions lie in a region very far from the perturbative regime, hence the fundamental 
QCD Lagrangian can not be responsible in this respect. Therefore, we need some non-perturbative approaches like 
QCD sum rules [l6| as one of the most powerful and applicable tools to hadron physics. Note that, the coupling 
constants, D*D S K, DtDK^, [lj], D Q DJK, D s0 DK D*Dp pU, D*Dtt (2(J ^LDDp (H, DDJ/iP [Hj], 
D*DJ/ijj [H], D*D*n [H D*D*J/i) |3, D S D*K, D* S DK DDuj [H, D*D*p [30|, and B s0 BK, B sl B*K 
[3l| have been investigated using different approaches. 

This paper is organized as follows. In section 2, we give the details of QCD sum rules for the considered coupling 
constants when both D(B) and K*(892) mesons in the final state are off-shell. The next section is devoted to the 
numerical analysis and discussion. 



II. QCD SUM RULES FOR THE COUPLING CONSTANTS 

In this section, we derive QCD sum rules for coupling constants. For this aim, we will evaluate the three-point 
correlation functions, 

n% B) (p\q) d\ d*y e*-* e^(0\T(vf(x) V D ^(y) ^^(O)) |0> (1) 

for D(B) off-shell, and 

n£;( 892 V,g) = t 2 J <fix d% e*'-* e^y(0\T(v DiB) (x) vf(v) ^^(O)) 10} (2) 

for if* (892) off-shell. Here T is the time ordering operator and q — p — p' is transferred momentum. Each meson 
interpolating field can be written in terms of the quark field operators as following form: 

r)f (x) = s(x)j v u(x) 
il D[B] {y) = u{y) lb c[b](y) 

v? [B:] (o) = s(o)iAm (3) 

The correlation functions are calculated in two different ways. In phenomenological or physical side, they are obtained 
in terms of hadronic parameters. In theoretical or QCD side, they are evaluated in terms of quark and gluon degrees 
of freedom by the help of the operator product expansion (OPE) in deep Euclidean region. The sum rules for the 
coupling constants are obtained equating the coefficient of a sufficient structure from both sides of the same correlation 
functions. To suppress contribution of the higher states and continuum, double Borel transformation with respect to 
the variables, p 2 and p' 2 is applied. 

First, let us focus on the calculation of the physical side of the first correlation function (Eq.([T])) for an off-shell 
D(B) meson. The physical part is obtained by saturating Eq. (JTJ) with the complete sets of appropriate D°, D* 
and A"* (892) states with the same quantum numbers as the corresponding interpolating currents. After performing 
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four-integrals over x and y, we obtain: 

(0\ v r\K*(p\e))(0\ V D ^\D(B)(q))(K%p\e)D(B)( q )\D:(B;^ 



il 2 - m2 D(B))(P 2 - m i t {B i ))(P' ~ m K') 

(4) 



where .... represents the contribution of the higher states and continuum. The matrix elements appearing in the 
above equation are defined in terms of hadronic parameters such as masses, leptonic decay constants and coupling 
constant, i.e., 

(0\ V f\K*(q,e))=m K »f K .e u 

(0\ V ^\D(BW)}=i ml ^[ D{B) 
m c {b) + m u 



*/ 



(D* s (B;)(p,e')\ V °° (BD \0) = m 2J;{fl;) / D:(B;) e 
(K*( qi e)D(B)(p')\D;(B;)(p,e')) = ig^K-^BK^^^pPc, (5) 

where 9d*d^K'(b*bk*) ^ s cou P nn g constant when D{B) is off-shell and e and e' are the polarization vectors associated 
with the K* and D*(B*), respectively. Using Eq. (|5|) in Eq. ^ and summing over polarization vectors via, 

q v qe 



e^eg = —g v e 



j , * 



K* 
PrjPfi 



S £ m = ~3w + ~2 . ( 6 ) 

" l D%{B* s ) 

the physical side of the correlation function for D(B) off-shell is obtained as: 

m 2 

nfD > , „ fD*(B*)fD(B)fK" m rrr rn D ,j B ,)m K * 

! V VP >P) — 9 d * dk *( B *bk')\1 > , 2 Za \7I5 Z5 Za T E PpPa + ---- [<) 



To calculate the coupling constant, we will choose the structure, from both sides of the correlation functions. 

From a similar way, we obtain the final expression of the physical side of the correlation function for an off-shell K* 
meson as: 



n /-D* (B*)fD(B)fK* ~^~7~~ TTljj* ( B *) m K* 
II* ( P ',P) = -g^DK.^BK^ ) , ,2 2 h 2 2 U 2 ^ VgP. + (8) 

' (P' -™ D{B ))(P - m D;(B:))(Q 2 - m K») 

Now, we concentrate our attention to calculate the QCD or theoretical side of the correlation functions in deep 
Euclidean space, where p 2 — > — oo and p' 2 — > — oo. For this aim, each correlation function, n^(p',p), where i stands 
for D(B) or K*, can be written in terms of perturbative and non-perturbative parts as: 

n;„ (p' , P ) = (n per + n nonper ) e a ^y' p p a , (9) 

where the perturbative part is defined in terms of double dispersion integral as: 

^[ds'fds ^* 2 ) 
4tt 2 J J (s - p 2 )(s' - p' 2 ) 

here, p(s,s',q 2 ) is called spectral density. In order to obtain the spectral density, we need to calculate the bare loop 
diagrams (a) and (b) in Fig. ([I]) for D(B) and K* off-shell, respectively. We calculate these diagrams in terms of 
the usual Feynman integral by the help of Cutkosky rules, i.e., by replacing the quark propagators with Dirac delta 
functions: q i^ m i — ► (—2Tri)5(q 2 — m 2 ). After some straightforward calculations, we obtain the spectral densities as 
following: 

p D {s, s',q 2 ) = . c - {2m 3 s q 2 + m u s (2to 2 - q 2 + s - s') - m 2 s m u (q 2 + s - s') + 2mjV + m 2 c [m s (-q 2 + s- a') 

A / ^5, 5 , Q ) 

- m u {-q 2 + s + s')] + m c [to 2 (-q 2 + s- s') - (q 2 + s - s') s' - to 2 (-q 2 + s + s')] 

q 2 + S -.s')+q 2 [-q 2 + s + s>)}}, (11) 



Ilper = j / ds' / ds - — - g — h subtraction terms, (10) 



(^,2 i 2 

to, to ' 
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FIG. 1. (a) and (b): Bare loop diagram for the D(B) and K* off-shell, respectively; (c) and (e): Diagrams corresponding to 
quark condensate for the D(B) off-shell; (d) and (f): Diagrams corresponding to quark condensate for the K* off-shell. 



pf (s, s', q 2 ) = vs/2 / >C g/ ^ {2m 3 c q 2 + m u s (2m 2 u - q 2 + s - s') - m 2 c m u (q 2 + s - s') + 2m? a s' + m 2 s (m c (-q 2 + s- s') 

(12) 



\ 3 / 2 (s,s',q 2 ) 

m u (-q 2 + s + s')) + m s [m 2 c {-q 2 + s - a') - (q 2 + s - s') s' - m 2 u (-q 2 + s + s')] 



m c [ml (q 2 +s- s') + q 2 (-q 2 +s + s')] } 



for the D*DK* (892) vertex associated with the off-shell D and K*(892) meson, respectively, and 

N n 



P a {s,s',q 2 ) 



A 3 / 2 (s,s',q 2 ) 

m u (-q 2 + s + s')] + m b [m 2 s (-q 2 + s- s') - (q 2 + s - s') s' - m 2 u (-q 2 + s + s')] 
™s {ml (q 2 + s- s') + q 2 [-q 2 + s + s')] } , 



P2 {s,s',q z ) = 



I „2\ 



Nr. 



{2m 3 s q 2 + m u s (2m 2 u - q 2 + s - s') - m 2 s m u {q 2 + s - s') + 2m\s' + ml [m s {-q 2 + s - s') 

(13) 

{2m 3 q 2 + m u s (2ml — <? 2 + s — s') — m\m u (q 2 + s — s') + 2m 3 s s' + m 2 (m& (— q 2 + s — s') 

(14) 



\y 2 (s,s',q 2 ) 

m u (~q 2 + s + s')) + m s [m 2 b (-q 2 + s - s') - (q 2 + s - s') s' - m 2 u (-q 2 +s + a')] 
m b [ml (q 2 + s- a') + q 2 (-q 2 + S + s')] } , 



for the B*BK* (892) vertex associated with the off-shell B and if* (892) meson, respectively. Here X(a, b, c) = a 2 + 
b 2 + c 2 — 2ac — 26c — 2ab and N c = 3 is the color number. 

To calculate the nonperturbative contributions in QCD side, we consider the quark condensate diagrams presented 
in (c), (d), (e) and (f) parts of Fig. (fTJ). It should be reminded that the heavy quark condensates contributions are 
suppressed by inverse of the heavy quark mass, so they can be safely neglected. Therefore, as nonperturbative part, 
we only encounter contributions coming from light quark condensates. Contributions of the diagrams (d), (e) and (f) 
in Fig. (fTJ) are zero since applying double Borel transformation with respect to the both variables p 2 and p' 2 will kill 
them because of appearing only one variable in the denominator in these cases. Hence, we calculate the diagram (c) 
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in Fig.([T]) for the off-shell D(B) meson. As a result, we obtain: 



D 



for the off-shell D meson and 



If 



nonper 



(p 2 — m 2 )(p f2 — mD 



(ss) 



(15) 



(16) 



for the off-shell B meson. 

Now, it is time to apply the double Borel transformations with respect to the p 2 (p 2 — > M 2 ) and p' 2 — > (p' 2 — > M' 2 ) 
to the physical as well as the QCD sides and equate the coefficient of the selected structure from two representations. 
Finally, we get the following sum rules for the corresponding coupling constant form factors: 



(q 2 - m 2 D ) 



-e m 2 e M ' 2 



!d*Jd!k> 



-m D «m K « 



4 TT 2 



ds 

(m c +m s ) 2 J (m s -\-m u ) 2 



da'p D (s,s',q 2 ) 



9[1 — (f D (s, s')) ]e«^e«' ! + (ss}e«^e«^ 



(17) 



(q 2 - m K * 2 ) 



D* D 

-e m' 2 eF 5 



!d»Jd!k> 



m c +m. 



-m D ,m K , 



1 



4 7T 



so 



ds 



(m c -\-m s ) 2 •* (m c -\-m u ) 



d S 'pf( S , S ',q 2 ) 



9[l-(ff(s,s')) }e^eTF 

for the off-shell D and if* (892) meson associated with the D* S DK* (892) vertex, respectively, and 

{q 2 - m%) 



-e m 2 e m' 



-m B >m K , 



4 TT 2 



2 - -a' m f ™ 2 

#[1 — (f B (s, s')) je^eM" + (ss)e"M 7 e"M T 



(m5+m s ) 2 (m s -t-m u ) 



dsV B (s,sV) 



(18) 



(19) 



(q 2 - m K , 2 ) 



Sb%SbSk* ^^m B: m K . 
6[l-{ff{a,s')f]e^e^ 



-e m' 2 e»" 



4 7T 



(is 

(mb+m s ) 2 J (rn b +m u ) 



ds'pf( S , S ',q 2 ) 



(20) 



for the off-shell B and if* (892) meson associated with the B*BK *(892) vertex, respectively. The integration regions 
in the perturbative part in Eqs. (|17 p - (|20p arc determined requiring that the arguments of the three S functions coming 
from Cutkosky rule vanish simultaneously. So, the physical regions in the a - a' plane are described by the following 
non-equalities: 



- 1< 



D = 2 s (to 2 - to 2 + s') + (m 2 c ~ to 2 - s)(-q 2 + s + s') < 1 

\ 1 / 2 (m 2 ,m 2 ,s)A 1 / 2 (s,s / ,q 2 ) 



1 </f M = 



2 s (-to 2 . + to 2 - s') + (to 2 - to 2 + s)(-q 2 + a + a') 



< 1, 



A 1 / 2 (to 2 , to 2 , s)A x / 2 (s, s',q 2 ) 
for the off-shell D and if* (892) meson associated with the D*DK* (892) vertex, respectively, and 



i</ B M') 



2 s (to 2 — to 2 + a') + (ml — ml — ■s)(-q 2 + a + s') 
A 1 / 2 (to 2 , to 2 , s)A 1 / 2 (s, s', g 2 ) 



< 1, 



-1</: 



r, ;x _ 2 s (-to 2 + to 2 - s') + (to 2 - to 2 + s)(-g 2 + a + s') 
2 1 ' ' X 1 / 2 (m 2 b ,m 2 ,s)X 1 / 2 (s,s',q 2 ) 



(21) 
(22) 

(23) 
(24) 
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for the off-shell B and K*(892) meson associated with the B*BK* (892) vertex, respectively. These physical regions 
are imposed by the limits on the integrals and step functions in the integrands of the sum rules. In order to subtract 
the contributions of the higher states and continuum, the quark- hadron duality assumption is used, i.e., it is assumed 
that, 

p higheratatea(^ ^ = p OPE ^ _ _ (25) 

Note that, the double Borel transformation used in calculations is defined as: 

B 1 - S. ( 1V"+" 1 1 r ~rn^/M 2 -m 2 JM' 2 1 / 2fi N 

(p 2 - m\) m (p' 2 - ml)" k ' T(m)T(n) ^ (M 2 ) m - 1 (Af 2 )"- 1 ' V ' 



III. NUMERICAL ANALYSIS 



Present section is devoted to the numerical analysis of the sum rules for the coupling constants. In further analysis, 
we use, m^.(892) = (0.89166 ± 0.00026) GeV, m D o = (1.8648 ± 0.00014) GeV, m D * = (2.1123 ± 0.0005) GeV, 
m B ± = (5.2792±0.0003) GeV , m B * = (5.4154±0.0014) GeV m c = 1.3 GeV, m b = 4.7 GeV, m s = 0.14 GeV|l, 
m u = 0, f K , = 225 MeT/JH, f D , = (272 ± 16°_ 20 ) MeV, f B . = (229±20i\ 6 ) MeV H, f B = (190 ± 13) MeV [H, 
f D = (206.7 ±8.9) MeV [37| and (ss) = -0.8(0.24 ± 0.01) 3 GeV 3 (H- 

The sum rules for the strong coupling constants contain also four auxiliary parameters, namely the Borel mass 
parameters, M 2 and M' 2 and the continuum thresholds, so and s' Q . Since these parameters are not physical quantities, 
our results should be independent of them. Therefore, we look for working regions at which the dependence of coupling 
constants on these auxiliary parameters are weak. The working regions for the Borel mass parameters M 2 and M' 2 
are determined requiring that both the contributions of the higher states and continuum are sufficiently suppressed 
and the contributions coming from higher dimensions are small. As a result, we obtain, 8 GeV 2 < M 2 < 25 GeV 2 
and 3 GeV 2 < M' 2 < 15 GeV 2 for D off-shell, and 4 GeV 2 < M 2 < 10 GeV 2 and 3 GeV 2 < M' 2 < 9 GeV 2 
for K* off-shell associated with the D*DK* (892) vertex. Similarly, the regions, 14 GeV 2 < M 2 < 30 GeV 2 and 
5 GeV 2 < M' 2 < 20 GeV 2 for B off-shell, and 5 GeV 2 < M 2 < 20 GeV 2 and 5 GeV 2 < M' 2 < 15 GeV 2 for 
K* off-shell are obtained for the B*BK* (892) vertex. The dependence of considered coupling constants on Borel 
parameters for different cases are shown in Figs.([2][5]) and ([71 1X0]) . From these figures, we see a good stability of the 
results with respect to the Borel mass parameters in the working regions. The continuum thresholds, sq and s' are 
not completely arbitrary but they are correlated to the energy of the first excited states with the same quantum 
numbers. Our numerical calculations lead to the following regions for the continuum thresholds in s and s' channels 
for different cases: (rnjy^B*) + 0.3) 2 < sq < (m D *( B ,) + 0.5) 2 in s channel for both off-shell cases and two vertexes, 
and {m D(B) + 0.3) 2 < s\ < (m D{B) + 0.7) 2 and + 0.3) 2 < s' < {m K * + 0.7) 2 for K* and D(B) off-shell cases, 

respectively in s' channel. Here, we should stress that the analysis of sum rules in our work is based on, so called the 
standard procedure in QCD sum rules, i.e., the continuum thresholds are independent of Borel mass parameters and 
q 2 . However, recently it is believed that the standard procedure does not render realistic errors and the continuum 
thresholds depend on Borel parameters and q 2 and this leads to some uncertainties (see for instance [Hj|). 

Now, using the working region for auxiliary parameters and other input parameters, we would like to discuss the 
behavior of the strong coupling constant form factors in terms of q 2 . In the case of off-shell D meson related to the 
D*DK* vertex, our numerical result is described well by the following mono-polar fit parametrization shown by the 
dashed line in Fig. ©: 

{D) -103-34 , } 

9d* s dk*W >- q 2 _ 28 57' ^ 7 > 

where Q 2 = —q 2 - The coupling constants are defined as the values of the form factors at Q 2 = —fn 2 neson (see also 
[l9j). where m meson is the mass of the on shell meson. Using Q 2 = —rr? D in Eq. (|2"T|) . the coupling constant for 
off-shell D is obtained as: g^, DK * = 3.23 GeV^ 1 . The result for an off-shell K* meson can be well fitted by the 
exponential parametrization presented by solid line in Fig. (|6]) , 

g%l K > (Q 2 ) = 4.44 e ™ - 0.70. (28) 

Using Q 2 — —m 2 K , in Eq. (f2"5|). the g^, DK , = 4.25 GeV -1 is obtained. Taking the average of two above obtained 
values, finally we get the value of the gu* dk* coupling constant as: 

g D , DK , = (3.74 ± 1.38) GeV' 1 . (29) 



7 



From figure (j6j) it is also clear that the form factor, g%, DK * is more stable comparing to g^-'DK-' with respect to 
the Q 2 . The similar observation has also obtained in [19( in analysis of the D*Dp vertex. In our case, the two form 
factors coincide at Q 2 = 0.1612 GeV 2 and have the value 3.64 GeV' 1 very close to the value obtained taking average 
of the coupling constants for two off-shell cases at Q 2 = —m 2 neson . 

Similarly, for B*BK* vertex, our result for B off-shell is better extrapolated by the mono-polar fit parametrization, 

-354.37 



SB) 



^BK 



Q 2 -98.14' 

presented by dashed line in Fig. (fTTj) and for K* off-shell case, the parametrization 



(30) 



gf»2 K , (Q 2 )= 3.02 e^r- 0.28, 



(31) 



shown by the solid line in Fig. (fTT|l . describes better the results in terms of Q 2 . Using Q 2 
coupling constant is obtained as g B * BK * = 2-78 GeV~ 
Eq. (|31~j) . Taking the average of these results, we get, 



-m% in Eq. (j30|), the 
Also, g B » BK , — 3.69 GeV' 1 is obtained at Q 2 = —vr? K , in 



g B , BK » = (3.24 ± 1.08) GeV~ 



(32) 



The errors in the results are due to the uncertainties in determination of the working regions for the auxiliary 
parameters as well as the errors in the input parameters. From the figure Fig. (|11[) . we also deduce that the heavier 
is the off-shell meson, the more stable is its coupling form factor in terms of Q 2 . From this figure, we also see that 
the two form factors related to the B*BK* vertex coincide at Q 2 = —0.7152 GeV 2 and have the value 3.58 GeV' 1 
also close to the value obtained taking average of the corresponding coupling constants for two off-shell cases at 
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FIG. 2. gE'DK- (Q 2 = 1 GeV 2 ) as a function of the Borel mass M 2 . The continuum thresholds, s = 6.83 GeV 2 , s' = 2.54 GeV 2 
and M' 2 = 5 GeV 2 have been used. 
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FIG. 3. gB*DK*(Q 2 ~ 1 GeV 2 ) as a function of the Borel mass M' 2 . The continuum thresholds, so = 6.83 GeV 2 , s' = 
2.54 GeV 2 and M 2 = 15 GeV 2 have been used. 
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FIG. 5. gE'DK»(Q 2 — 1 GeV 2 ) as a function of the Borel mass M' 2 . The continuum thresholds, so = 6.83 GeV 2 , s' = 
6.57 GeV 2 and M 2 — 5 GeV 2 have been used. 
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FIG. 6. Qd* dk* as a function of Q 
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•BK'iQ 2 — 1 GeV 2 ) as a function of the Borel mass M 2 . The continuum thresholds, so = 34.99 GeV 2 , s' 
2.54 GeV 2 and M' 2 = 10 GeV 2 have been used. 
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FIG. 8. <?s*Bif*(Q 2 = 1 GeV 2 ) as a function of the Borel mass M' 2 . The continuum thresholds, so = 34.99 GeV 2 , s' 
2.54 GeV 2 and M 2 = 20 GeV 2 have been used. 
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FIG. 9. 9b* bk* iQ 2 = 1 GeV 2 ) as a function of the Borel mass M 2 . The continuum thresholds, s = 34.99 GeV 2 , s' 



35.75 GeV 2 and M' 



GeV have been used. 
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FIG. 10. gs'BK'iQ 2 = 1 GeV 2 ) as a function of the Borel mass M' 2 . The continuum thresholds, so = 34.99 GeV 2 , 
s' = 35.75 GeV 2 and M 2 — 8 GeV 2 have been used. 




FIG. 11. gB'BK" as a function of Q 2 . 



